Observation of a Pressure-Induced Topological Quantum Phase Transition in BiTel 
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We report the observation of a pressure-induced topological quantum phase transition in the polar 
semiconductor BiTel using X-ray powder diffraction and infrared spectroscopy. The X-ray data 
confirm that BiTel remains in its ambient-pressure structure up to 8 GPa. The lattice parameter 
ratio c/a shows a minimum between 2.0-2.9 GPa, indicating an enhanced c-axis bonding through p z 
band crossing as expected during the transition. Over the same pressure range, the infrared spectra 
reveal a maximum Drude spectral weight, reflecting the closing and reopening of the semiconducting 
gap through band inversion. Both of these features are characteristics of a topological quantum phase 
transition, and are consistent with a recent theoretical proposal. 

PACS numbers: 64.70.Tg, 71.70.Ej, 78.20.-e, 61.05.cp 
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Topological insulators are a class of materials with an 
insulating bulk and protected metallic edge or surface 
states [l[ . This intriguing state of matter was theoreti- 
cally predicted and experimentally observed 
in a variety of 2D and 3D systems. In the search for 
topological insulators, a guiding principle is to identify 
materials having band inversion at special time-reversal- 
invariant momentum points of the Brillouin zone, or to 
induce such a band inversion in suitable systems 0, Q . 
This necessary condition was convincingly demonstrated 
inBiTKSx^Se^ [HHland (Bij.Jn^Seg 0,113 by 
doping, and in bismuth [18J by substrate-induced strain. 
The effects of doping and strain in these experiments are 
to tune the spin-orbit interaction and/or the lattice pa- 
rameters, driving the system through a topological quan- 
tum phase transition (TQPT). Such effects can be real- 
ized directly by pressure, a method that does not in- 
volve the potential defects and inhomogeneity of dop- 
ing and offers more tunability than substrate-induced 
strain. This has indeed been proposed as a means to 
look for a topological insulating state in, e.g., ternary 
Heusler compounds [l9j |. A^hoO^ (A = Y or rare-earth 



elements) [20|, Ge 2 Sb 2 Te 5 [2l|, and BiTel [22j. Despite 



its importance to the field, there has been no experimen- 
tal indication for a pressure-tuned topological insulator. 

The prediction of pressure-induced TQPT in BiTel is 
particularly interesting 22]. BiTel is a layered polar 
semiconductor with the trigonal space group P3ml and 
a hexagonal unit cell at ambient pressure [23[ . The triple 
Te-Bi-I layers stack along the crystallographic c axis, 
coupled by the Van der Waals interaction. The broken 
inversion symmetry in the presence of strong spin-orbit 
coupling gives rise to a giant Rashba-type spin splitting, 



the largest observed so far in any system [24j. Although 
BiTel is a narrow-gap semiconductor with an ambient- 
pressure band gap of 



-0.38 eV 24-2 



it is typically 

self-doped due to non-stoichiometry and thus an n-type 
semiconductor 



271 ] . Remarkably, Bahramy et al. 



predicted that hydrostatic pressure effectively tunes the 
crystal-held splitting and spin-orbit interaction in BiTel, 
turning it into a topological insulator. Such a TQPT is 
accompanied by a band inversion near the A point of the 
Brillouin zone; the dispersion of the bulk spin-split bands 
becomes almost linear at the critical pressure P c . 

In this Letter, we report experimental observation of 
characteristics of such a pressure- induced TQPT in BiTel 
using X-ray powder diffraction and infrared spectroscopy. 
X-ray powder diffraction shows that BiTel remains in the 
hexagonal structure up to 8 GPa but with the lattice pa- 
rameter ratio c/a passing through a minimum between 
2.0 and 2.9 GPa, indicating an enhanced c-axis bonding 
through p z band crossing. Infrared spectroscopy over the 
same pressure range reveals a maximum in the free car- 
rier spectral weight as would be expected when the semi- 
conducting band gap closes and reopens through band 
inversion. Both of these results provide evidence for a 
pressure-induced TQPT in BiTel at a critcal pressure P c 
between 2.0 and 2.9 GPa [22|. Our experiments demon- 
strate pressure as an efficient way to induce TQPTs and 
the use of infrared spectroscopy and X-ray diffraction for 
their investigation. 

Single crystals of BiTel were grown by the Bridgman 
method. The ambient-pressure carrier density is esti- 
mated to be 2.7xl0 19 cm -3 [28[. Samples used in dif- 
ferent experiments were mechanically cleaved from one 
single crystal. High-pressure infrared spectroscopy and 
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FIG. 1. (color online), (a) X-ray powder diffraction pattens 
of BiTel from 0.9 to 9.7 GPa. The data above 0.9 GPa are 
vertically shifted for clarity. The filled circles in (b), (c), and 
(d) are the lattice parameters c and a, the unit-cell volume, 
and the c/a ratio in the hexagonal phase at different pres- 
sures, extracted from our XRD data. The open circles are 
the ambient-pressure data taken from Ref [23l |. The solid line 
in (c) is a fit to a Birch-Murnaghan equation of state. The 
error bars in (d) come from the uncertainties in c and a. (e) 
A diagram illustrates the band inversion across the TQPT 
at P c . Magenta denotes the Bi-6p z orbital; cyan denotes the 
Te-5p z and l-5p z orbitals. 



angle-dispersive X-ray powder diffraction were performed 
at U2A and XI 7C beamlines at the National Synchrotron 
Light Source (Brookhavcn National Laboratory). All 
data were collected at room temperature. Pressure was 
monitored by laser-excited ruby fluorescence. The sep- 
aration of the two fluorescence lines indicated that the 
pressure was quasi-hydrostatic. In the X-ray diffrac- 
tion (XRD) experiment, a sample was ground into fine 
powder and loaded into a diamond anvil cell (DAC). A 
4:1 methanol-ethanol mixture was used as the pressure- 
transmitting medium. Two-dimensional diffraction rings 
were collected for pressures up to 27.2 GPa, with the in- 
cident monochromatic X-ray wavelength set to 0.4066 A. 
Integrating the diffraction rings yielded XRD patterns as 
a function of the diffraction angle 29 [29j |. 

The experimental configuration for the infrared mea- 
surements is illustrated in the inset of Fig. [3]b. The sam- 
ple was pressed against a diamond anvil and surrounded 
by a pressure-transmitting medium (KBr) elsewhere. In- 



frared microspectroscopy was performed on an FT-IR 
spectrometer coupled to a microscope (30j . Near-normal 
incidence reflectance and transmittance were measured 
for photon energies between 0.08-1.00 eV and pressures 
up to 25.2 GPa. Four samples were measured, all giv- 
ing consistent results. The data with the most complete 
pressure dependence are presented here. 

The XRD patterns in Fig. show a structural 
phase transition at ~8 GPa. Ambient-pressure BiTel 
has a hexagonal lattice with lattice parameters a = 
4.339 A and c = 6.854 A [H]. Increasing pressure consis- 
tently shifts each Bragg peak to a larger 29 angle, indi- 
cating the reduction of the d-spacing. Above 7.9 GPa, a 
great number of new peaks emerge, suggesting a different 
structure. We analyzed the XRD patterns [3l[ and dis- 
cuss here the results for the hexagonal phase. The anal- 
ysis yields the lattice parameters (Fig. [TL) and the unit- 
cell volume (Fig. [Tfc). We include the ambient- pressure 
values of these parameters taken from Ref [23| (omitted 
in our experiment), shown as open circles in Fig. [T|3,c; 
they are consistent with the pressure dependence of our 
data. A least-square fit of the unit-cell volume to a Birch- 
Murnaghan equation of state [HH yields the bulk modulus 
Bo = 23.1 GPa and its pressure derivative B' = 7.2. 

The ratio c/a shown in Fig. [T]i has a clear minimum 
between 2.0-2.9 GPa, strongly indicating a change in the 
chemical bonding across the TQPT. In a topologically- 
trivial insulator, a TQPT occurs when the atomic orbital 
energy ordering near the Fermi level reverses and the con- 
duction (valence) band inverts. The diagram in Fig. \Tjp 
illustrates this process. In ambient-pressure BiTel, the 
bottom conduction bands are dominated by the Bi-6p z 
orbital while the top valence bands are dominated by the 
Te-5p z (and l-5p z ) orbital 22, 33]; that is, within the low- 



energy bands, the bonding is almost ionic in nature for 
P < P c . After the pressure-induced band inversion hap- 
pens (P > P c ), one expects the Te-5p z orbital contributes 
to the bottom conduction bands and the Hi-6p z orbital 
dominates the top valence bands; the bonding becomes 
more covalent. At the transition (P = P c ), additional 
metallic bonding occurs involving the p z orbitals of Bi 
and Te, whose charge fluctuation would make the c-axis 
more compressible. This is observed clearly as a min- 
imum in the ratio c/a, which indicates a TQPT at P c 
between 2.0-2.9 GPa. The corresponding volume con- 
traction is 7-9%, similar to the theoretically predicted 
11% [22j |. Thus, we have the first piece of strong evi- 
dence of a TQPT. 

Next, infrared spectroscopy reveals more detailed in- 
formation on the electronic structure in this pressure 
range and indicates a TQPT as well. In the hexagonal 
phase, the reflectance (Fig.[2ji) shows a plasma edge char- 
acteristic of the free carrier response in the system. In 
0-2.20 GPa, the plasma edge gradually blue shifts upon 
pressure increase. Above 2.20 GPa, it plunges to a lower 
photon energy and keeps red shifting when the pressure 
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FIG. 2. (color online). Infrared reflectance (a) and transmit- 
tance (b) of BiTel for pressure in 0.17-7.51 GPa. The arrows 
indicate the direction of increasing pressure. Absorption from 
diamond corrupts the data between 0.22-0.33 eV (patched as 
dashed lines) and causes the dip at ~0.46 eV in reflectance. 
The inset in (b) shows the measurement configuration. The 
arrowed lines indicate the beam paths for the reflected (R) 
and transmitted (T) signals. 



is further increased. The transmittance (Fig. ^jp) also 
shows dramatic pressure dependence. Between 0.08- 
0.20 eV, it changes slowly below 2.20 GPa but jumps 
to a higher level at 2.45 GPa, indicating a reduction of 
absorption. Both the reflectance and the transmittance 
contain fringes; they arise because of the Fabry-Perot ef- 
fects in the optically-flat sample and in the KBr. The 
thinner sample gives rise to the broad fringes and the 
thicker KBr gives rise to the fine ones [3(J. When pres- 
sure is increased above 2.20 GPa, the contrast of the 
broad fringes in reflectance improves significantly. This 
suggests a sudden reduction of absorption right above 
2.20 GPa at the photon energies where the broad fringes 
reside, corroborating the aforementioned sudden increase 
in transmittance. Evidently a transition happens in the 
2.20-2.45 GPa range, coinciding with the pressure where 
the minimum c/a occurs. We note that the turning of 
the plasma edge at ~2 GPa was consistently observed in 
the four samples we measured, although the quantitative 
behavior of the infrared data is sample dependent on ac- 
count of differing scattering rates in the samples. For the 
discussion that follows and to simplify the analysis, we 
remove the fringes by Fourier transforming the data and 



FIG. 3. (color online), (a) Reflectance at low pressure and 
low photon energy, after fringe removal. The solid lines are 
fits to the Drude model. The inset illustrates the conduc- 
tion band dispersion e(k) and the chemical potential as the 
pressure approaches P c . The slope of e(k) at the chemical po- 
tential gives the Fermi velocity vf- The arrows indicate the 
direction of increasing pressure, (b-d) Pressure dependence 
of the fitting parameters for the Drude model, i.e. the plasma 
frequency uj p , the electronic scattering rate 7, and the dielec- 
tric constant eoo due to high-energy incoherent excitations. 
The solid lines in (b) are linear fits to the parts before and 
after the maximum. 

filtering the fringe signatures. 

We adopt the following approach to analyze the in- 
frared spectra after fringe removal. The sample and 
the pressure-transmitting medium are considered as 
a double-layer slab sandwiched between the diamond 
anvils. Incoherently adding intensities of multiple inter- 
nal reflections in the double layers gives the reflectance 
and transmittance [3fj| . The equations involve the sample 
thickness, which was obtained by estimating its ambient- 
pressure value and using the pressure dependence of c 
determined from XRD 30] . 

Below 0.2 eV, the free carrier absorption dominates. 
The Drude model adequately describes the free carrier 
response, with a dielectric function 



Here is the dielectric constant contributed from high- 
energy incoherent excitations, uj p the plasma frequency, 
and 7 the electronic scattering rate. A least-square fit to 
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Eq. (SI) in [3(j was performed on reflectance, shown as 
solid lines in Fig. [3^,. The fitting parameters are shown in 
Fig. [3b,c,d. The plasma frequency at the lowest pressure 
is 5.59xl0 14 s _1 , close to the ambient-pressure value ob- 
tained from an infrared measurement performed on the 
same sample [26|. Using the effective mass m* = 0.1 9mo 
(mo is the electron mass) (28|, we estimate the free car- 
rier density n — 1.9 x 10 19 cm~ 3 from cu p = yj Anne 2 / m* , 
reasonably consistent with that estimated from quantum 
oscillation measurements [2^. Compared to an earlier 
ambient-pressure measurement 26j, the scattering rate 



7 is bigger in the sample studied here, and varies among 
all the samples measured. The monotonic increase in 
upon pressurization accords with that in reflectance at 
photon energies greater than 0.3 eV. 

The peak at 2.20 GPa in the pressure-dependent lu p 
shown in Fig. [5b is a second piece of evidence for the 
TQPT. Across the TQPT in an insulator, the band gap 
closes and reopens. Near the band touching point, the 
dispersion of the conduction and valence bands becomes 
nearly linear. A simplified picture is illustrated in the in- 
set of Fig. [3Ji. If electrons occupy the conduction band, 
the change of the conduction-band dispersion causes a 
change of the Fermi velocity vf- This manifests as a 
change of the Drude spectral weight cjp/8 originating 
from the carriers lying close to the chemical potential. 
As illustrated in the inset of Fig. [3^, vf increases as 
P —> P c (and decreases above P c when the dispersion 
again becomes parabolic, not shown). The maximum 
vf lies in the transition region where the band disper- 
sion approaches linear behavior, yielding the maximum 
uj p around the transition because w 2 oc v 2 F [Hj]. More- 
over, vf is expected to be linear in pressure to the lowest 
order 3J], resulting in the linear dependence of <jJ p (P) 
before and after the transition, shown as solid lines in 
Fig. \Sjp. Thus the observed maximum in the Drude spec- 
tral weight clearly indicates the closing and reopening of 
the semiconducting gap in our naturally doped sample 
during a TQPT; the critical pressure is consistent with 
the range given by the X-ray powder diffraction. 

We confirm our analysis by the optical conductivity. 
Using Eqs. (SI) and (S2) in [3(J, we iteratively solve 
for the sample refractive index and extinction coefficient. 
These are used to calculate the optical conductivity; its 
real part is shown in Fig. 0] The data may be quan- 
titatively inaccurate because of the analysis invoked to 
account for the complicated measurement configuration, 
but they show the intraband and interband transitions 
reported previously in the 0.1-0.6 eV range [HI, (2(| and 
hence should be qualitatively correct. The dc conduc- 
tivity o o , roughly estimated by extrapolating the Drude 
conductivity to zero photon energy, peaks at 2.20 GPa, 
consistent with the maximum in gq = uj^/Att^i calcu- 
lated using the ui p and 7 from the Drude fit (inset in 
Fig. 0}. The sudden drop of the Drude conductivity 
above 2.20 GPa results in the fringe contrast improve- 
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FIG. 4. (color online). The real part of the optical con- 
ductivity calculated from reflectance and transmittance. The 
arrow indicates the direction of increasing pressure. The data 
between 0.22-0.33 are corrupted by the diamond absorption 
(patched as dashed lines). The inset shows the dc conductiv- 
ity (To calculated from ui p and 7 shown in Fig. [3b, c. 



ment in reflectance and the increase in transmittance 
shown in Fig. [5] The band gap closing and reopening 
associated with the TQPT is not discernible as interband 
features. This could be accounted for by the self-doping 
in our sample: the band gap change is severely masked 
by the Pauli blocking effects because free carriers occupy 
the conduction band. Measurements on BiTel samples 
with intrinsic semiconducting properties, such as that re- 
ported in Ref 27|, will provide more useful information 
of the band gap. 

In summary, we obtained independent and consistent 
experimental evidence for a pressure-induced topological 
quantum phase transition in BiTel using X-ray powder 
diffraction and infrared spectroscopy. Across the transi- 
tion: (i) The lattice parameter ratio c/a shows a min- 
imum; (ii) The Drude spectral weight peaks. Both of 
these features are characteristics of the topological quan- 
tum phase transition, and are consistent with a recent 
theoretical proposal. Our work motivates the exploration 
of the novel physics associated with this transition 35 1 
using other experimental probes. 
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Supplemental Material 



Infrared microspectroscopy 
in a diamond anvil cell 



A diamond anvil cell (DAC) with type Ha diamond 
anvils was used in the infrared microspectroscopy. In- 
frared spectra were collected on a Bruker Vertex 80v FT- 
IR spectrometer coupled to a Hyperion 2000 microscope, 
using a Globar source and an MCT detector. The mea- 
surement configuration is shown in Fig. IS 1 1 The sam- 
ple was pressed against the top diamond anvil and sur- 
rounded by a pressure-transmitting medium (KBr) else- 
where. The sample covered approximately half of the 
300 (im diameter sample chamber in a stainless steel 
gasket. BiTel is soft; when compressed in a DAC, it 
forms a shiny flat surface with the top diamond anvil 
to make specular reflections. At a given pressure, we 
first obtained the transmittance: the infrared beam was 
focused on the KBr next to the sample to collect a ref- 
erence transmission spectrum T ri and then on the sam- 
ple to collect the sample transmission spectrum T s ; the 
ratio T s /T r gave the sample transmittance. We next ob- 
tained the reflectance: the infrared beam was first fo- 
cused on the sample to collect the sample reflection spec- 
trum R s , and then on the top diamond- air interface to 
collect a reference reflection spectrum R r , yielding a ra- 
tio Ri = R s /R r . When all the pressure-dependent mea- 
surements were finished, we opened the DAC and took 
the top diamond anvil, which had been in contact with 
the sample. After cleaning its culet, we measured the 
reflection spectrum from the top diamond-air interface 
(R r ) and from the culet-air interface through the anvil 
(i? c ) to yield another ratio R2 = R c /R r - The sample 
reflectance at the sample-diamond interface was calcu- 
lated as R1/R2 ■ Rd — R s /Rc ■ Rd, where Rd is the known 
single-bounce reflectance of diamond at the diamond-air 
interface. Such a measurement scheme minimized possi- 
ble systematic errors and the effects due to the absorption 
from diamond. 

The spectra contain broad fringes due to the Fabry- 
Perot effects in the sample and fine fringes due to the 
same effects in the pressure-transmitting medium KBr. 
Because of our measurement scheme, both broad and fine 
fringes are present in the reflectance, but only fine fringes 
dominate in the transmittance. For reflectance, the sam- 
ple spectrum R s contains the Fabry-Perot effects from 
both the sample and KBr, but the spectra R r and R c 
do not have fringes. The ratio R s /R c therefore preserves 
both the broad and the fine fringes. Because the sample 
is absorbing, and because the beam has to traverse the 
sample twice to make the Fabry-Perot effects in KBr ob- 
servable in the measured reflectance, the corresponding 
fringe contrast is weak. Therefore these fine fringes are 
seen as small oscillations superimposed on a background 



of the much broader sample fringes. For transmittance, 
the reference spectrum T r was taken in the region next to 
the sample, where KBr is thicker than that underneath 
the sample by an amount equal to the sample thickness. 
This thickness mismatch gives rise to the dominant fine 
fringes in the transmittance. 



2. Reflectance and transmittance 
of a double-layer slab 

The measured reflectance and transmittance can be 
analyzed by considering the sample (denoted as s) and 
the pressure-transmitting medium (denoted as m) as a 
double-layer slab, shown in Fig.EU This slab is bounded 
by the diamond anvils (denoted as d) on both sides. Since 
the fringes were removed from the spectra to simplify the 
analysis, one can add intensities rather than amplitudes 
from the multiple internal reflections in s and to. Because 
the reflectance and transmittance of a single-layer slab is 
straightforward to obtain, we solve the double-layer slab 
problem following a two-step method 0. (i) Solve the 
single-layer-slab problem, where the slab is to, bounded 
by s and the bottom d. This gives the reflectance R s . m d 
and transmittance T s . m d = 1 — R s -md of a composite 
medium (denoted as md) consisting of m and the bot- 
tom d, interfaced with s. (ii) Solve the single-layer-slab 
problem, where the slab is s, bounded by the top d and 
the composite medium md. The reflectance and trans- 
mittance for the double-layer slab are 
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FIG. SI. (a) Measurement configuration for infrared mi- 
crospectroscopy in a diamond anvil cell. The arrowed lines 
indicate the beam paths for the reflected (R) and transmit- 
ted (T) signals, (b) Close-up of the part in the square in 
(a), d: diamond anvil, a: sample, m: pressure-transmitting 
medium, md: composite medium consisting of m and d. 
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In Eqs. (3), (JS1, and (JS3|, = [(n-n d ) 2 + K 2 ]/[(n + 
n d f + k 2 ], R sm = [(n - n m ) 2 + K 2 ]/[(n + n m f + k 2 ], 
and R m d = (n rn — n d ) 2 /(n m + n^) 2 are the single-bounce 
reflectance at the s-d, s-m, and m-d interfaces, respec- 
tively. They involve the to-be-solved sample refractive 
index n and extinction coefficient k, and the known re- 
fractive index of diamond Q and KBr n m (sj. is 
the wavenumber (defined as 1 /wavelength with the unit 
cm -1 ), and x is the sample thickness. The determination 
of x is discussed in the next section. 

3. Sample thickness 

The thickness of the sample in the DAC can be de- 
termined from the fringes in the reflectance spectra. 
We Fourier transformed the raw reflectance spectrum at 
0.17 GPa and filtered out the fringe signature due to 
KBr. This makes the broad fringes due to the sample 
easier too identify. The broad fringes have a period of 
Av « 300 cm -1 . Taking the refractive index of the sam- 
ple calculated from the Drude fit at this pressure in this 
photon energy range, n « 4.1, we determined the sample 
thickness x — \j2nAv w 4.1 /im. Note that the sample 
refractive index n is not a known quantity. We initially 
fitted the reflectance to the Drude model using the single- 
bounce reflectance formula for the diamond-sample inter- 
face, R s d = [{n - n d ) 2 + K 2 ]/[(n + n d ) 2 + k 2 ]. The refrac- 
tive index n derived from the fit was used to calculate the 
sample thickness. The Drude fit was then refined using 



Eq. (|Sip . We checked the resulting n and repeated these 
procedures, until a self-consistent solution was found. To 
determine the sample thickness at other pressures, one 
can use the same procedures described above, or use the 
pressure dependence of the lattice parameter c obtained 
from X-ray powder diffraction. The latter approach was 
adopted in our analysis because it is simpler and proba- 
bly more accurate compared to the former. 

To confirm the validity of the thickness estimation 
based on fringes, we applied this method to determine 
the thickness of the pressure-transmitting medium KBr. 
We obtained the fine-fringe period from the raw trans- 
mittance spectrum at 0.17 GPa, Av « 83 cm -1 . The 
refractive index of KBr is about 1.53 in the photon en- 
ergy range where the fine fringes reside. Therefore the 
substrate thickness is roughly 1/(2 x 1.53 x 83 cm -1 ) = 
39 /xm. Note that we measured the thickness of the 
stainless-steel gasket to be ~41 nm before we compressed 
it between the diamond anvils. This is reasonably close 
to the estimated KBr thickness. 
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